Fermi Surface Nesting and the Origin of the Charge Density Wave in NbSe2 
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We use highly accurate density functional calculations to study the band structure and Fermi 
surfaces of NbSe2. We calculate the real part of the non-interacting susceptibility, Kxo(q), which is 
the relevant quantity for a charge density wave (CDW) instability and the imaginary part,SJxo(Q)i 
which directly shows Fermi surface (FS) nesting. We show that there are very weak peaks in 5Rxo(q) 
near the CDW wave vector, but that no such peaks are visible in Sfxo(q)5 definitively eliminating FS 
nesting as a factor in CDW formation. Because the peak in 5Rxo(q) is broad and shallow, it is unlikely 
to be the direct cause of the CDW instability. We briefly address the possibility that electron- 
electron interactions (local field effects) produce additional structure in the total (renormalized) 
susceptibility, and we discuss the role of electron-ion matrix elements. 
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In 1964 V.L. Ginzburg proposed excitonic supercon- 
ductivity in quasi-2D structuresi composed of metal 
layers sandwiched between insulating layers. By that 
time, only two layered materials were known to be 
superconducting,^-^ PdTe2 and NbSe2. In the four 
decades since then, NbSe2 and isostructural selenides 
have been intensively investigated (close to 1.5 thousand 
papers published to date). However, the main interest in 
this compound has shifted from the fact that it is a lay- 
ered material (hundreds of layered superconductors are 
now known), to the existence of a nearly-commensurate 
charge density wave (CDW>i instability and it's possi- 
ble interplay with the superconductivity that sets in at 
a lower temperature. We mainly forego discussion of the 
interesting issues surrounding the superconducting state, 
its origin, and its relationship to the CDW state, and in- 
stead concentrate on the mechanism behind the CDW 
transition itself. 

The first electronic structure calculations for NbSe2 
were presented by Mattheiss in 1973^. Using a non-self- 
consistent potential he was able to produce a band struc- 
ture with basic features in reasonable agreement with 
more recent self-consistent calculations^'^-^, but which 
showed only two bands crossing the Fermi energy (it is 
now known that there are three), and underestimated the 
energy depth of a saddle point at ~ ^ TK. Fermi surfaces 
based on this band structure led to early suggestions that 
the CDW transition was driven by nesting^, an assump- 
tion that has carried through to the present time. The 
nearness of the saddle point to the Fermi energy (Ei?) 
led Rice and Scotli^ to argue that CDW formation was 
driven, not by Fermi surface (FS) nesting in the con- 
ventional sense, but rather by saddle points lying within 
ksTcDW oi Fjp and separated by the CDW wavevec- 
tor, QcDW = (5,0,0). A significant amount of effort 
has been spent on resolving the 'controversy' between 
the nesting and saddle point theories for NbSe2 and for 
related CDW compounds such as 2H-TaSe2, lT-TaSe2, 
TaS2 and others, but no specific feature that would give 
rise to an instability at Qcdw has been convincingly 
isolated. As early as 1978, Doran et aJifi^ used a sim- 
ple one-band tight-binding model to show that the sus- 
ceptibility has no sharp peaks of electronic origin and 



made a rudimentary calculation of the electron-ion ma- 
trix elements that suggested strong coupling at Qcdw- 
Nonetheless, in the more modern era of self-consistent 
calculations, the FS nesting and saddle point hypothe- 
ses continue to be debated and FS nesting is cited as a 
contributing factor in CDW formation in almost every 
paper that addresses the issue. In this paper we report a 
quantitative analysis of the non-interacting susceptibility, 
Xo{^.,'^ = 0), based on first principles density functional 
theory (DFT) calculations, and show that not only is the 
structure in momentum space extremely weak, but that 
FS nesting does not play any role at all and can be defini- 
tively ruled out as a cause for CDW formation. Instead, 
we argue, the instability must be due to the electron- 
ion interaction, or possibly from local field effects. The 
structure of the paper is as follows: In Section I, we lay 
out the basic mathematical framework for dielectric cal- 
culations, and discuss the various approximations we em- 
ploy. In Section II, we present and discuss the calculated 
band structure and Fermi surfaces and show that there 
is good agreement with experiment. In Section III, we 
show the real and imaginary parts of the non-interacting 
susceptibility, and discuss the origins of prominent fea- 
tures, emphasizing the lack of a strong peak at the CDW 
wave vector and the irrelevance of FS nesting. 



I. GENERAL THEORY 

A material is unstable toward the formation of a CDW 
at wave vector Q when the response function (inverse 
electronic dielectric function) fails to satisfy the stability 
conditiorki^, 

e-i(Q,o) = i + yc(Q)x(Q,o)<i (1) 

where Vc{q) = 47re^/q^ is the Coulomb potential, and 
the susceptibility x(q, 0) is normally negative. A dielec- 
tric constant within the range < e(Q,<jj) < 1 therefore 
signals a CDW instability and, in conjunction with Eq. 
n indicates that the latter range is accessible only from 
the lower boundary, i.e when x(Q,a;) diverges. In the 
simplest approximation, x(Qi is calculated by ignor- 
ing the exchange part of the electron-electron interaction 
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in favor of Hartree-like terms only. This is known as the 
random phase approximation (RPA)^'^ and relates x to 
the non- interacting susceptiblity, xo^ 

xrpa{(i, t^) = xo(q, w)/[i - Vb((j)xo(q, w)], (2) 

where xo is 



|(Zk|e''i''|Z',k + q)| 



(3) 

Here f{s) is the Fermi function and the matrix ele- 
ments (Zk|e*'i'"|/', k + q) are between Bloch functions, 
l^k) = u;k(r)e''"'. If these do not have much structure 
as a function of q, one can set them uniformly to unity, 
yielding the following expression for the imaginary part 
of Xo: 



lim 3xo(q,t^)/^ 



n 



(2-)^ 



|V/k X V;',k+q| 



(4) 



where the line Lk is the intersection of the two surfaces 
defined by ejk = and £i'k+q = (the Fermi energy is 
set to zero). Eq. 0] correctly suggests that a Fermi sur- 
face nesting, that is, a Fermi surface topology such that 
the length of the line Lk is particularly large at some 
Q, should lead to a peak in lim^^^o Q'XoCQ, '^)/'^- This 
is correct, but misses several important points. First, as 
Eq. ^ indicates, a strong maximum also requires small 
and/or nearly parallel vjk and v;/ k+Q along ik- Second, 
it is the real, and not the imaginary part of the static sus- 
ceptibility that defines a CDW instabihty (at w = 0, Sxo 
vanishes). In contrast to 3xoj which indeed depends only 
on electronic characteristics near the Fermi surface, Kxo 
collects information from an energy window of the order 
of the band width (which follows from Eq. |3| taking into 
account that J dE/E diverges at large arguments). 
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Finally, one should keep in mind that a divergence of 
the electronic susceptibility signals an instability of the 
system with respect to spontaneous formation of a CDW 
even if the nuclei are clamped to their high-symmetry po- 
sitions. Of course, if released, the nuclei would shift as 
well. However, divergence of x is a substantially more 
severe criterion than the soft mode condition, which is 
that the frequency of a particular phonon mode at some 
wave vector becomes zero. This latter condition defines 
an instability toward the freezing in of a particular pe- 
riodic pattern of ionic displacement, which is the only 
CDW that is actually observed in real systems. 

Phonon frequencies are related to the dielectric func- 
tion via the so-called Pick-Cohen-Martin formulaic. Here 
we give this formula in a simplified form for an elemental 
solid: 



Da^piq) = ^„/3(q)-^a/?(0), 

^aMq) = E lq+G|'^c(q + G)(q + G)„ X (5) 
e-l(q + G, q -f G')(q + G%Fc(q + G'), 



where Daf3{ci) is the dynamical matrix and a, P are 
Cartesian indices. It is clear that an eigenvalue of this 
matrix can soften and eventually become zero even if the 
macroscopic dielectric functioniS*i&: 



e(q)-l/[e-l(q+G,q + G')]c 



(6) 



does not diverge. Here and above, the dielectric function 
e"i(q + G, q G') differs fro m e^jp^ in that it includes 
the so-called local field eff'ects: Umklapp processes and 
the exchange-correlation interaction between electrons. 
Formally, it is written (neglecting for simplicity the band 
indices) as: 



e-i(q + G,q+G')= ^gg' + ^ ^C7(q + G)xo(q + G, q + G^) ® 

Gi 

[Sa,G' -E{'5GiG.^c(q + Gi)-/xc(q+Gi,q+G2)}xo(q + G2,q+G')]-^ (7) 
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(in these expressions, matrix inversion is with respect to 
the reciprocal lattice vector indices), where Ixc{^,^') = 
(5^£^xc/^p(r)'5p(r')is the exchange correlation kernel 
[note that in the LDA /xc(i" - r') = (5(r - r')/xc(i"), 
and thus /xc(q G, q + G') - Ixc{G - G')] . 

Therefore, several separate questions must be posed: 
1) Does the geometry of the Fermi surface provide nesting 
in the sense of a strong maximum of Eq^? 2) Is there a 
peak in the real part of xoj as defined in Eq|31at the wave 



vector corresponding to the observed CDW and if so, is 
it related to any nesting-derived structure in 3x(qiO)? 
3) Is there a divergence in e(q) as defined by Eq|H| ? 4) 
If the previous answer is no, is there a soft mode (zero 
eigenvalue) in Eq. 5? 

Below we present numerical density functional calcu- 
lations that provide answers to the first two questions 
and supply information allowing for reasonable conjec- 
ture about answers to the last two. The substance of our 



3 



findings is that there is only weak nesting in this system, 
and at the "wrong" wave vector. While Kxo does have 
a broad maximum at the CDW wave vector, it is likely 
too small to account for the CDW instability. Lacking 
full linear response calculations, we cannot exclude with 
certainty an electronic instability in e(q) (Eq. EJ, due to 
local field effects (exchange-correlation and Umklapp) be- 
yond plain RPA, but it is more likely that the instability 
appears only after electron-ion interactions are explicitly 
taken into account (Eq. 5). Experiment further supports 
our contention as only the phonon corresponding to the 
actual CDW is observed to soften, whereas, if e(q) were 
to diverge at some q, all phonons would soften at this 
vector. 

To illustrate the last point we note that xrpa itself 
can have a peak, but obviously cannot diverge (Eq. [2J|. 
Furthermore, including the exchange-correlation interac- 
tion between electrons within DFT, but without Umk- 
lapp processes, gives 

= [l + {Vfc(g)-/xc(g)}Xo(q,0)] 

This expression cannot diverge either. Finally, it can 
be shown that, generally speaking, the Umklapp pro- 
cesses tend to lower the electronic dielectric constant^^. 
Since e(q, w) can become unstable only by passing from 
a negative quantity to a positive one, there can be no 
CDW formation at any value of q. Non-local correla- 
tion effects neglected by DFT may be instrumental in 
creating an instability beyond Eq. |S1 but it is physically 
more likely that the electron-ion matrix elements, which 
can naturally have strong q-dependence, account for the 
CDW instability through their effect on the dynamical 
matrix. In fact, this is consistent with some early conjec- 
tures based on simplified tight-binding calculationsiitii. 
In these works, a strong phonon softening was observed 
at the correct wave vector, but only when electron-ion in- 
teractions were accounted for. The authors of both works 
conclude, correctly, that this softening is not simply re- 
lated to FS geometry, but stay short of excluding weak 
FS nesting as an essential factor. As we will show, no rel- 
evant nesting exists in NbSe2, and, therefore, it cannot 
play any role in the CDW instability. 

II. ELECTRONIC STRUCTURE 




FIG. 1: The crystal structure of NbSe2, showing the two dis- 
tinct layers that make up the unit cell. Each Nb ion sits in 
a hexagonal mirror plane with three Se neighbors above and 
below, forming a network of distorted edge-sharing octahedra. 

of NbSe2. It is comprised of two formula units arranged 
in layers, each layer containing a hexagonal Nb plane 
sandwiched between two shifted hexagonal Se planes. 
We used the experimentally measured lattice constantsSS, 
(a=3.44 A and c=12.55 A), which agrees well with other 
reports^LSSiS^ and we relaxed the Se ions to their lowest 
energy positions along the c axis (z=0. 1183c). This po- 
sition is in good agreement with all experimental reports 
(z = 0.116 - 0.118c)^V^i^^ except Ref. HI which gives 
z = 0.125c. Our relaxation resulted in a Se-Se interpla- 
nar distance of 3.57 A which can be compared with the 
Nb-Se intralayer distance of 2.59 A. Based on these sep- 
arations, we expect an anisotropic band structure with 
rather strong interplanar coupling. To calculate the non- 
interacting susceptibility of Eq. [SJ we used a direct sum- 
mation over a very fine mesh of 73x73x7 ^ 40,000 k 
points in the BZ, and a Fermi temperature smearing of 
T= 10 meV. 



Our calculations were carried out using a full-potential, 
augmented plane wave plus local orbital (APW-l-lo) 
schemeiS,. The exchange-correlation potential was ap- 
proximated by a local density approximation (LDA) 
parameterization^^ and the muffin-tin radii were set to 
2.3 (Nb) and 2.2 (Se) with a value of 7.0 for RKmax. We 
found small but important differences between calcula- 
tions with and without the spin-orbit interaction. Ex- 
cept where explicitly noted, we include spin-orbit cou- 
pling using a second variational method for all calcula- 
tions. The structure shown in Fig. ^displays the unit cell 



A. Band Structure 

Our band structure (Fig. |2) exhibits some differences 
from previous full potential calculations jSi^ giving rise to 
differences in the resulting Fermi surfaces. The discrep- 
ancies are most likely due to a large difference in Se po- 
sitions: z=0.134c in Refs. 0@ compared to our relaxed 
height of z=0.118c. We were not able to identify an ex- 
perimental paper reporting such high Se positions. 

We observe two classes of bands near the Fermi energy 
- the first is the fully antibonding Se Pz band, which. 
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FIG. 2: (color online) The band structure of NbSe2. Top: A 
closeup of the bands around the Fermi energy, neglecting spin- 
orbit coupling. The highlighted bands are of Se character 
and cross the Nb bands without hybridization. Bottom: The 
full band structure, including the spin-orbit interaction which 
relieves the Nb band degeneracy in the fez = tt/c plane. The 
highlighted bands are of combined x'^ ~ /xy character 



by parity, does not hybridize with the z— even d-orbitals 
of Nb {dx2-y2, dxy, and dz^-i) in the kz—0 plane, and 
does not hybridize with the z— odd Nb d orbitals in the 
kz = tt/c plane. In Fig. 12^7 a close-up of the band 
structure without spin-orbit interactions is shown. The 
Pz band is highlighted with filled circles and can be seen 
crossing and re-crossing the Nb band in the T — K — M 
plane without mixing, despite a nearly perfect energetic 
degeneracy. In the A — L — H plane, the range of the 
figure no longer includes the pz band which has dropped 
nearly 2 eV in energy. This rather dramatic energy low- 
ering can be traced to bonding-antibonding interactions 
between unit cells. The width of the Se p band complex 
(-6.0 eV to 0.5 eV) is determined by the energy difference 
between the bonding and anti-bonding configurations of 
the four pz orbitals along the c-axis within the (two layer) 
unit cell and between the cells themselves. In the kz—0 
plane, antibonding pz orbitals in adjacent unit cells are 
antibonding with each other, pushing the energy of the 
single fully anti-bonding band up to the Fermi level. In 
the kz = tt/c plane, the relationship between adjacent- 
cell Pz orbitals becomes bonding, the band is lowered, 
and the Fermi crossing is removed. Fig. O shows the an- 
tibonding configuration of the Pz orbitals at the F and 
A points, with arrows indicating the crucial interactions 
between cells. 

The second class of bands consists of one Nb d band. 



which, like the Pz orbitals, is split in the kz=Q plane, 
by bonding-antibonding interactions both within and be- 
tween unit cells. In Fig. |21 spheres representing the 
Nb dz^-i orbital (or equally the other z-even orbitals, 
— and xy) are shown in fully bonding and fully 
antibonding configurations at the F point. By contrast, 
in the kz = tt/c plane [A point shown), the bonding- 
antibonding energy difference arising from orbital orien- 
tations within the unit cell is counteracted by interactions 
between cells and the band therefore remains degenerate 
(see Fig. |3i)- In Fig- the full band structure, includ- 
ing spin-orbit, is shown. The spin-orbit coupling allows 
for p^-Nb mixing along V — K — M and removes the strict 
degeneracy of the Nb bands at k^ — it jc. It is tempting 
to identify the latter bands as simply Nb dz2_i orbitals, 
but is misleading: symmetry arguments, confirmed by 
the numerical decomposition of the LAPW bands, pre- 
dict that this band should have atomic d^^^i character 
along the FA line, but d.j.2_y2 -\- dxy character along the 
KH line and mixed character elsewhere. This can be 
seen clearly in Fig|2l3 the where the dx2_y2 -\- dxy bands 
are highlighted and are clearly dominant near the Fermi 
level along some symmetry directions. On can, of course, 
enforce mapping of this band structure onto a one-band 
dz2 — 1 TB model24, but this does not reflect the charac- 
ter decomposition of the LAPW bands nor does it mean 
that this band is physically derived from actual Nb dz2-i 
orbitals. 



B. Fermi Surfaces 

In Fig. 0]our three calculated Fermi surfaces are dis- 
played in an expanded (four total zones) Brillouin zone 
scheme. These surfaces are significantly different from 
early calculations , ^i^^i^^ most notably in the existence 
of a third band crossing Ep- On the other hand, they 
are quite similar to the sketches of Ref. 0, with the 
difference that we designate all three surfaces as "hole- 
like" because the filled states are external to the drawn 
surfaces, whereas in Ref. the third surface is called 
"electron-like" , presumably because the band is more 
than half-filled. Below we compare each FS to exper- 
imental data from angle-resolved photoemission spec- 
troscopy (ARPES) and de Haas van Alphen (dHvA) 
quantum oscillation measurements. 

Interlayer coupling breaks the degeneracy of the high- 
est unfilled Nb d band and gives rise to two concentric 
cylindrical Fermi surfaces around the F — ^ and K — L 
lines of the BZ. The surfaces derived from the higher 
(anti-bonding) energy band (Fig. 2J;) are almost dis- 
persionlcss along the kz direction, while those derived 
from the lower (bonding) energy band (Fig. ^p) are 
strongly warped. Neither of these have been seen with 
dHvA, which requires very pure samples to detect large 
orbits, but they have been consistently seen by ARPES 
measurements^iSiSLSSiSS which are in very good general 
agreement with calculation. We will use the measure- 
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FIG. 3: (color online) A schematic showing the origin of the 
splitting of the antibonding pz band and the Nb (1^2 _i band 
at two different kz positions. The Nb ions are omitted from 
the structural graphic for clarity and the (1^2 _i orbitals have 
been idealized as spheres. Left side: At the F-pt, Nb dz2_i 
orbitals are either fully bonding or fully antibonding, and the 
antibonding orbitals are antibonding even between cells 
(blue arrows). Right side: At the A-pt, there is no overall 
energy difference between bonding and antibonding Nb dz2_i 
configurations and the antibonding pz orbitals are now bond- 
ing between cells. Plus signs indicate bonding and minus 
signs antibonding - the long black arrows connect the top- 
most orbital to the one that would sit above it if the cells 
were repeated endlessly. 



ments of Rossnagel et a&, who investigated directions 
both parallel and perpendicular to the a-b plane, for spe- 
cific comparison with our surfaces. The observed Fermi 
surfaces derived from the higher Nb band (Fig. ^p) show 
less kz dispersion than in our calculations and conse- 
quently the pockets around K — L remain rounded along 
their full length rather than becoming sharply triangu- 
lar as ours do. The outer pockets of the lower Nb band 
(Fig. are not warped, consistent with our results, but 
the K — L centered triangular pockets again appear to 
be more rounded. Overall there is substantial agreement 
between ARPES measurements and theory, particularly 
in terms of the observed placement and kp oi the Nb- 
derived FS pockets - the essential elements of interest for 
Fermi surface nesting. 

The small pancake surface in Fig. is derived from a 
Se 4p band, and has an extremely small width along the 
kz direction. Very early dHvA measurements'^'^ saw this 
surface, but misattributed it by comparing to a band 
structure:^^ that didn't show the Se band crossing the 
Fermi energy. Later measurements by Onuki et aJ2i con- 
firm the existence of these surfaces, though the authors 
were still apparently unaware of the calculated Se-derived 
Fermi surface and therefore report a discrepancy with 
theory. The most recent dHvA experimental clearly de- 




FIG. 4: The three Fermi surfaces of NbSe2 in an extended 
BZ scheme with the conventional BZ indicated as a hexagon 
in the center of each figure, a) The small Se-derived (band 1) 
pancake surface around the F -point, b) The bonding (band 
2) Nb-derived Fermi surface with strong warping along kz- 
c) The nearly 2D anti-bonding (band 3) Nb-derived Fermi 
surface. The V — K — M plane is at the top and bottom of 
each figure and the A — L — H plane cuts through center. 



tect and correctly identify the smallest surface, obtain- 
ing good quantitative agreement with the earlier works. 
ARPES measurements see a weak signal of the smallest 
surfac o^i^^ with an in-plane area that agrees well with 
dHvA measurements, but with greatly exaggerated per- 
pendicular dimension. Our calculated Se Fermi surface 
is larger both in- and out-of-plane with respect to dHvA 
measurements, but shows a smaller perpendicular width 
than ARPES, which is not expected to be as accurate 
in this direction. For comparison with available experi- 
mental data, our calculated the dHvA frequencies for the 
Se-derived pocket as well as all the other Fermi surfaces 
is shown in Tabled 

Note that our calculated cross-sectional areas are no- 
ticeably different from previous work by Corcoran et al. 
This can be partially ascribed to the different Se po- 
sition in their calculations. However, there were also 
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TABLE I: The calculated dHvA frequencies of NbSe2. Bands 
two and three have two separate surfaces, designated in the 
Table by the axis that runs through their centers F — A or 
K — L. The two extremal orbits for each surface in a (001) 
field occur in the T — K — M and A — L — H planes, indicated 
by r and A respectively in the Table. For the Se pancake 
surface, the frequency was also calculated for an in-plane field 
direction (120). 



probably too few fc-points to ensure an accurate area 
calculation. Because of the simple shape of the Fermi 
surfaces, one can calculate their volumes from the ex- 
tremal cros-section, assuming a full cylindrical symme- 
try for the Se pocket (which yields V ~ Sxz^/T^Sxy), 
and using a trapezoidal integration for the open pockets 
(F « [Si + S2]t^/c). Applying this procedure to Corco- 
ran's numbers we obtain a violation of the Luttinger the- 
orem of the order of more than 1%, while our methodp22. 
gives the correct total number of electrons with an error 
of less than 0.01%. 

While ARPES data on the Nb Fermi surfaces agrees 
well with our calculations, the de Haas- van Alphen area 
for the Se pocket (the only one measured so far) differs 
by nearly 4000 T for the largest orbit and 500 T for the 
smallest one. However, due to the large effective mass 
for this pocket, a relatively minor charge transfer of only 
0.02 electrons from this band to the Nb bands is needed 
to bring the cross-sections into agreement. The geometry 
of the large Nb cylinders would not be changed in any 
important way with such a transfer. Therefore, we can 
safely use the calculated band structure to evaluate the 
nesting effects and susceptibility in NbSe2. 

The small pancake FS, although it carries little DOS 
and is probably not important for superconductivity or 
CDW formation, plays an important role in transport. 
We have calculated the Fermi velocities and plasma fre- 
quencies for both in-plane and perpendicular directions. 
A decomposition by band, as shown in Table ^J, reveals 
that while the pancake (band 1) contributes barely 5% 
of the total DOS, it is carrying more than three quarters 
of the total current in the c direction, and is responsible 



FS sheet 




DOS 


^px ^pz 


VFx 


VFz 






(states/Ry) (eV) (eV) (10 


* cm" 


{10^ cm 


Band 1 F - 


A 


3.8 


0.40 2.16 


0.10 


0.53 


Rand 2 F — 


A 


9.9 


1.63 0.78 


0.23 


0.11 


Band 2 K - 


-L 


19.7 


1.65 0.86 


0.18 


0.09 


Band 3 F - 


A 


11.4 


1.60 0.13 


0.22 


0.02 


Band 3 K ~ 


L 


29.5 


1.85 0.26 


0.16 


0.02 


total 




74.2 


3.37 2.46 


0.19 


0.14 



TABLE II: The density of states, plasma frequencies and 
Fermi velocities for the 5 sheets of the Fermi surface. 



for the relatively 3D character of the resistivity in this 
compound. If, as discussed above, we adjust the relative 
position of the Se pancake and the other bands, ojp^ for 
this band would drop by less than 50%, so that it would 
still contribute about as much as all other bands together 
to the transport across the planes. This, has important 
ramifications for superconductivity: any experiment re- 
lated to tunneling (Josephson contacts, Andreev reflec- 
tion etc..) in the direction perpendicular to the plane will 
probe mostly the Se pancake band which carries a very 
small part of the Cooper pairs and may actually have 
a reduced superconducting gap (reference to ARPES). 
Furthermore, if there is any difference between the su- 
perconducting properties of different Nb-derived sheets, 
this may also manifest itself in a nontrivial way in such 
tunneling, since the third band has nearly zero kz disper- 
sion and and its contribution to the c-transport is nearly 
neghgible (See Table HIJ. 

We now discuss the scenarios previously suggested 
as candidates for nesting at Qcdw and compare them 
quantitatively with our theoretical Fermi surfaces. The 
six symmetry-equivalent saddle points along the F — 
K directions (not shown) are separated by Qspi = 
(^, ^,0) or by Qsp2 = (5, 0,0), neither of which is near 
QcDW- Moreover, the calculated saddle point is approx- 
imately 150 meV beneath the Fermi energy, and although 
ARPES measurements^ show that band renormalization 
reduces the depth to only 50 meV, this is still too great an 
energy compared to the fcs7bDH'=3meV required for the 
saddle point theory of Rice and Scott^, as has been noted 
by other A There have been suggestions that self-nesting 
between the parallel flat edges of the central hexagonal 
Fermi surfaces around T — A could be responsible for 
the CDW instabilitjii. Our calculations show that only 
the smaller (second band) surface is oriented correctly to 
produce self-nesting along T—M; the larger surface would 
self-nest along T—K. The spanning vector between the 
faces of the inner surface is (0.41, 0, 0) which is too large 
in comparison to Qcdw (that of the outer surface is, of 
course, even larger) . A similar conclusion was reached 
in Refs. H and (2^ Furthermore, due to substantial kz 
dispersion in the second pocket, good nesting does not 
occur along the full length of the cylinder, substantially 
reducing Lk in Eq. 0] and thereby diminishing the con- 
tribution to X- Rossnagel et a&, noting this dispersion 



along kz, suggested that the broad peak calculatediiiii 
in XoiQcDW^O) is due to imperfect (weak) nesting be- 
tween these warped cyhndrical surfaces themselves (self- 
nesting) and between neighboring cylinders. Our calcu- 
lation of x(q, 0) shows that such a scenario is not man- 
ifested for NbSe2 (see Section III). It has been recently 
postulated-^ that the K — L centered cylinders, which are 
smaller than those around T — A, might nest at approxi- 
mately the right wave vector. The cross-sections of these 
cylinders, however, are not circles but are rounded trian- 
gles with flat edges oriented at 120° with respect to one 
another, making such nesting highly unlikely. The only 
apparent strong nesting possibility, according to our cal- 
culations occurs between the flat triangular edges of the 
K — L cylinders and the parallel flat edges of the central 
hexagonal pocket, Q = (i, |,0). In the next section, we 
will show that it is indeed this nesting that produces the 
highest peaks in the imaginary part of the susceptiblity. 
We reiterate here that FS nesting causes peaks in the 
imaginary part of xo: while the real part of xqj which is 
the essential quantity for CDW formation, draws from a 
much larger energy range and may have features entirely 
unrelated to FS topology. 



III. NON-INTERACTING SUSCEPTIBILITY 

The calculated real and imaginary parts of the non- 
interacting susceptibility are shown in Fig. |5l The matrix 
elements of Eq. Ohave been set to du> so that no inter- 
band contributions are taken into account. The effects 
of these terms turn out to be quite weak as we will show 
later in the text. The clear peaks in 3xo indicate nesting 
at q = (i, i, 0) and by symmetry, at q = (|, |, 0). These 
correspond to translating the central pocket from F to 
K. This results in partial nesting between the hexagonal 
and triangular pockets, as well as between two triangu- 
lar pockets with different orientations. Since the pockets 
can slide around inside one another, the peak is broad 
and, in fact, is not truly a single peak, but a compos- 
ite peak composed of three separate peaks near q = _fC, 
and corresponding to separate nestings of the three edges 
of each triangle. Assuming that the triangular surfaces 
are, in actuality, more rounded as indicated by ARPES, 
we expect these three peaks to further "smear" into one 
central peak. There is no indication of any other nesting 
and, most significantly, there is no peak whatsoever at 
q = QcDW- This situation differs from that of isostruc- 
tural NbTe2 where a strong peak in 3x(q, 0) does exist^ 
at QcDW, though we note that such a peak cannot be 
taken as evidence that nesting contributes to CDW for- 
mation without a calculation of 5Rx(Qi 0)- Unless the true 
fermiology of NbSe2 differs signficantly from our calcula- 
tions (and comparison with experiment strongly suggests 
that it does not), there is no evidence of a FS nesting 
contribution to CDW formation at all. 

The peaks in the imaginary part of the susceptibil- 
ity are always present in the real part because Ep is 




QcDW 



FIG. 5: The non-interacting susceptibility of NbSe2. a) The 
imaginary part exhibits FS nesting driven peaks at Q = 
(^,^,0). The plane at the bottom is a guide for the eye 
and corresponds to the lowest value of ^Xo- b) The real part 
has very weak peaks at Qcdw that come from energy inter- 
vals away from Ef. (see text). The plane corresponds to the 
lowest value of SRxo , which is also the density of states at the 
Fermi level. 



part of the sampled energy range, but these can be ei- 
ther emphasized or overshadowed by contributions from 
other energies. In the case of NbSe2, the peaks at K 
are no longer even visible and broad maxima appear in- 
stead along the T — M line. These do occur approxi- 
mately at QcDWj but are unlikely to be strong enough 
to produce a CDW instability. A two-dimensional plot 
of 5Rxo along F — M is shown in Fig. El with the posi- 
tion of the observed CDW wave vector indicated. The 
highest point along this curve is 33.5 Ryd~^ above the 
lowest point, 3fJxo(0)=N(Ei?)=70.3 Ryd"^, and occurs at 
q=(0.31,0,0) which is, in fact, rather near Qcdw- Fig. 
includes a band-by-band decomposition of the contri- 
butions to 5Rxo(q) along F — M showing almost no con- 
tribution at all from the Se-derived band and a rather 
peaked contribution from the antibonding Nb band, but 
at a wave vector quite different from Qcdw- Finally, 
we have included a plot of 5)?Xo(q) with both intra- and 
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FIG. 6: Kxo(q) along the T-M line of the Brillouin zone. The 
total is shown in shaded grey with the contributions from each 
individual band delineated by different lines. The height of 
band 1 has been multiplied by two to make it visible above the 
X-axis. The topmost line includes non-diagonal (interband) 
contributions. 



interband matrix elements set to unity in Eq. This ex- 
aggerates the contributions of interband transitions, but 
gives a rough idea of their effect on the susceptibihty. We 
find an overall shift upward of the entire spectrum but 
no particular sharpening of the peaks. Along T — M , our 
susceptibility somewhat resembles that of Ref. 0, but 
the suggestion that the broad peak is due to self-nesting 
of the central pocket is clearly unjustified as no such peak 
appears in Sxo- 

To isolate the true source of the (weak) maxima, we 
calculate a function, 



(9) 



from which the susceptibility can be recovered through 

'■^ F{x,y) 



dx / dy- 

-oo Jep X-y 



(10) 



By adjusting the limits of the integrals in Eq. ^1 the 
energy interval that gives rise to any given feature in xo 
can be isolated. For a very small interval near the Fermi 
energy, all the features of Sxo are reproduced, including 



the dominant feature at K . However, for energy intervals 
away from the Fermi energy, there is little or no contri- 
bution to xo {K) , with the result that the real part of xo 
has a minimum there when all energies are taken into 
account. Much of the contribution to SRxo comes from 
eigenvalues in the energy range |ek-i-q — ek| > 0.02 Ryd. 
This reinforces in a rather dramatic way our earlier as- 
sertion that one should not expect a CDW instability to 
stem from FS nesting. 

IV. CONCLUSIONS 

We have calculated the real and imaginary parts of the 
non-interacting susceptibility as a function of wave vec- 
tor throughout the Brillouin zone. The real part of xo 
exhibits very weak peaks near Qcuw that are unlikely to 
be strong enough alone to cause a CDW instability. The 
main contributions to these peaks come from an energy 
range not near 'Eip^ indicating that FS nesting is irrele- 
vant. A direct calculation of 9xo(q) confirms that there 
is absolutely no FS nesting at Q,cdw, although peaks 
due to nesting along q = (i, i,0) are evident. It can 
therefore be definitively stated that Fermi surface nesting 
contributes nothing to CDW instability in NbSe2. This is 
indicative of the more general, but often overlooked fact 
that Fermi surface nesting has an immediate effect only 
on the imaginary part of susceptibility, relevant, for in- 
stance, for the spin excitation spectrum, but not to CDW 
formation. The latter is defined by the structure of the 
real part of susceptibility, which is not directly affected 
by the Fermi surface nesting. More precisely, the effect 
of nesting on KxoCQ)'^) occurs in a very small part of 
the total relevant energy range (which is of the order of 
electronvolts rather than Kelvins), and if the electronic 
structure changes qualitatively within this range for a 
given system, which is nearly always the case, then nest- 
ing cannot be considered a primary mechanism of CDW 
formation. 
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